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Progress  has  been  made  In  the  program  to  use  solar  seismology  as  a  tool  for  the  understanding 
and  potentially  the  forcasting  of  solar  variability.  Changes  In  the  velocity  of  sound  profile  and 
nonadlabatic  effects  in  the  internal  equilibrium  condition  of  the  convection  zone  over  a  six  year 
period  have  been  inferred  from  observed  changes  in  the  frequencies  of  low-order  acoustic  moses 
and  Intermediate  degree  f-modes.  These  changes  may  be  associated  with  the  eleven  year  solar 
cycle.  The  observed  spectrum  of  gravity-modes  has  been  used  to  infer  a  discontinuity  in  the  mean 
molecular  weight  In  the  Sun’s  core  and  to  place  an  upper  limit  on  the  density  of  WIMPS  (Weakly 
Interacting  Massive  Particles)  proposed  to  explain  the  solar  neutrino  paradox.  The  interred  upper 
limit.  If  confirmed  in  future  work,  removes  the  WIMPS  proposal  as  a  viable  option.  Evidence  for 
gravity-mode  coupling  has  been  extended  from  4  multiplets  to  20.  This  evidence  for  mode 
coupling  leads  to  interesting  predictions  concerning  the  solar  neutrino  flux.  New  differential  radius 
observations  were  obtained  in  1988  and  1989.  Lrc,  wj- 

*7  • 


irwmTTT,. 


Sun;  Solar  Seigmology;  Global  Oscillations  '  1  j 


17.  SICUNTY  CLASSIFICATION 

UK&tfSPiED 


09  AUTIUCT 

UNCLASSIFIED 


AFOSR  Contract  No.  AF0SR-87-01 05 


APPLICATION  OF  SOLAR  SEISMOLOGY  TO 
STUDY  OF  THE  CONVECTION  ZONE 


Henry  A.  Hill,  Principal  Investigator 
Department  of  Physics 
University  of  Arizona 
Tucson,  Arizona  85721 
(602)  621-6782 


30  April  1990 

Final  Technical  Report  for  January  1,  1987  to  December 


Prepared  for 

THE  AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 


Ms.  Kathleen  L.  Wetherell,  Contracting  Officer 
Dr.  Henry  R.  Radoski,  Program  Manager 
AFOSR/NR 

Bolling  Air  Force  Base 
Washington,  DC  20332-6448 


,  1989 


TABLE  OF  CONTENTS 


Page 


Summary . 3 

1.  Research  Progress  From  1  Jan.  through  31  Dec.  1989 

Under  AFOSR  Contract  No.  87-0105 . 4 

2.  Publications . 10 

3.  Advanced  Degrees  Awarded . 12 


4.  1988-89  Oral  Presentations  presented  at  Conferences  and  Seminars  .  .13 

Attachment:  Report  Documentation  Page  (DD  Form  1473,  83  APR)  .  .  .14 


2 


Summary 


Progress  has  been  made  in  the  program  to  use  solar  seismology  as  a  tool 
for  the  understanding  and  potentially  the  forecasting  of  solar  variability. 
Changes  in  the  velocity  of  sound  profile  and  ncnadiabatic  effects  in  the 
internal  equilibrium  condition  of  the  convection  zone  over  a  six  year  period 
have  been  inferred  from  observed  changes  in  the  frequencies  of  low-order 
acoustic  modes  and  intermediate  degree  f-modes.  These  changes  may  be 
associated  with  the  eleven  year  solar  cycle.  The  observed  spectrum  of 
gravity-modes  has  been  used  to  infer  a  discontinuity  in  the  mean  molecular 
weight  in  the  Sun's  core  and  to  place  an  upper  limit  on  the  density  of  WIMPS 
(Weakly  Interacting  Massive  Particles)  proposed  to  explain  the  solar 
neutrino  paradox.  The  inferred  upper  limit,  if  confirmed  in  future  work, 
removes  the  WIMPS  proposal  as  a  viable  option.  Evidence  for  gravity-mode 
coupling  has  been  extended  from  4  multiplets  to  20.  This  evidence  for  mode 
coupling  leads  to  interesting  predictions  concerning  the  solar  neutrino 
paradox  and  temporal  variations  in  the  solar  neutrino  flux.  New 
differential  radius  observations  were  obtained  in  1988  and  1989. 
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1.  Research  Progress  From  1  January  1988  to  31  December  1989 
Under  Contract  No.  AFOSR-87-01 05 

The  research  objective  is  to  study  solar  variability  through 
observation  of  certain  types  of  variability  and  to  examine  properties  of  the 
solar  convection  zone  relevant  to  understanding  the  driving  mechanisms  of 
active  regions.  The  ability  to  anticipate  certain  types  of  variability  will 
likely  depend  on  an  understanding  of  these  mechanisms. 

The  observational  approach  is  two-fold:  (1)  direct  observation  of 
certain  types  of  solar  variability  —  changes  in  the  solar  diameter,  solar 
oblateness,  and  limb-darkening  function  --  that  are  accessible  through  solar 
diameter  observations,  and  (2)  the  study  of  the  solar  convection  zone  with 
solar  seismology.  In  regards  to  the  latter,  the  discovery  of  global 
oscillations  of  the  Sun  has  made  it  possible  to  study  the  spatial  and 
temporal  behavior  of  the  equilibrium  conditions  of  the  convection  zone. 

The  research  has  been  involved  with  the  development  of  techniques  for 
making  observations,  data  reduction  and  data  interpretation  and  the 
acquisition  of  data  that  will  permit  the  study  of  the  ei genf  requency 
spectrum  of  the  Sun. 

* 

Work  was  conducted  at  SCLERA  on  the  analysis  of  total  solar  irradiance 
observations.  An  objective  has  been  to  better  understand  the  properties  of 
gravity  mode  signals  in  total  irradiance.  It  was  found  that  a  significant 
fraction  of  the  gravity  mode  power  in  total  irradiance  observations  appears 
in  sidebands  of  the  classified  gravity  mode  frequencies.  These  sidebands 
appear  to  be  harmonics  of  the  rotational  frequencies  of  the  nonuniform  solar 
surface.  This  discovery  will  probably  be  quite  valuable  in  addressing  such 
questions  as:  (1)  Why  are  total  irradiance  gravity  mode  signals  so  small 
relative  to  those  obtained  in  differential  radius  observations  at  SCLERA? 
and  (2)  Why  do  various  groups  report  highly  discrepant  results  for  the 
internal  rotation  rate  of  the  Sun?  These  findings  were  presented  at  the 
XXVII  Plenary  COSPAR  Meeting,  July  18-29,  1988  in  Helsinki,  Finland  (Kroll, 
Hill  and  Chen  1988). 

A  second  paper  that  appeared  in  the  COSPAR  proceedings  contains  the 
most  recent  SCLERA  results  on  changes  over  a  six-year  period  in  the  low- 
order,  low-degree  acoustic  and  f-mode  spectrum  of  the  Sun  (Oglesby,  Hill  and 
Rabaey  1988).  These  particular  sets  of  modes  are  most  sensitive  to  the 
Sun's  convection  zone  and  therefore,  the  observed  changes  in  frequencies  of 
oscillation  are  very  likely  due  to  changes  in  the  convection  zone. 


# 

SCLERA  is  an  acronym  for  the  Santa  Catalina  Laboratory  for  Experimental 
Relativity  by  Astrometry,  a  research  and  educational  program  operating  from 
both  the  Solar  Studies  Division  and  the  Department  of  Physics  in  the  Faculty 
of  Science  at  the  University  of  Arizona. 
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Temporal  shifts  of  the  referenced  low-order,  low-degree  p-mode  and 
intermediate-degree  f-mode  eigenf requencies  motivated  a  sensitivity  analysis 
of  the  ei genf r equenci es  in  preparation  for  a  formal  inversion  (Rosenwald, 
Hill  and  Czarnowski  1988).  The  most  interesting  result  from  this  analysis 
is  that  the  frequencies  for  the  f-modes  possess  a  relatively  high 
sensitivity  to  the  entropy  eigenfunction  6s  in  the  convection  zone  in 
comparison  to  that  for  the  low-order  acoustic  modes.  Estimates  of  implied 
changes  in  the  internal  structure  of  the  Sun  have  been  obtained  using  the 
results  of  the  sensitivity  analysis.  It  was  inf°^red  that  Ac/c  -  0.01%  in 
the  upper  half  of  the  convection  zone,  where  c  is  the  sound  speed  (Oglesby 
and  Hill  1988)  . 

The  work  on  the  extension  of  the  range  of  gravity-mode  classification 
was  submitted  for  publication  in  1988  and  has  been  published  (Hill  and  Gu 
1990).  The  results  reported  in  this  paper  were  successfully  tested  by  using 
the  SCLERA  1978  diameter  observations  (Hill  1988),  the  SCLERA  1983 
differential  radius  observations  (Yi  and  Hill  1988),  and  the  1986-87 
differential  radius  observations  (Cornuelle  and  Hill  1989).  The 
classifications  were  also  tested  in  a  fourth  independent  analysis.  It  has 
been  predicted  by  Provost  and  Berthomieu  (1988),  that  the  distribution  of 
observed  amplitudes  of  the  temperature  eigenfunction  of  low-degree  gravity 
modes  should  exhibit  a  minimum  for  ei genf requenci es  corresponding  to 
frequencies  of  divergence-free  f-mode  like  oscillations  at  the  surface.  The 
amplitudes  of  the  gravity  modes  classified  by  Hill  and  Gu  (1990)  were 
examined  for  minima  at  the  predicted  frequencies.  The  plot  of  power  density 
versus  frequency  for  the  l  -  3  modes  has  a  well-defined  minimum  at  v  =  1 48 
pHz  while  the  plots  for  X,  =  2  and  4  modes  have  changes  in  slopes  at  v  =  115 
uHz  and  =  1 80  uHz,  respectively  (cf.  Hill  et  al .  1988).  The  presence  of 
these  3  observationally  obtained  frequency  features  near  the  theoretically 
predicted  locations  of  Provost  and  Berthomieu  give  support  to  the 
interpretation  that  the  observed  signals  are  due  to  oscillations  in  the 
solar  atmosphere  and  that  the  l  classifications  of  the  signals  are  also 
correct . 

The  asymptotic  properties  jf  the  low-degree  gravity  modes  classified  by 
Hill  and  Gu  (1990)  were  studied  in  the  framework  of  asymptotic  theory 
predictions  (Gu,  Hill  and  Rosenwald  1990).  One  of  the  more  interesting 
findings  of  this  analysis  places  a  significant  upper  limit  on  the  particle 
density  of  WIMPS  which  have  been  proposed  to  explain  the  solar  neutrino 
paradox.  The  inferred  upper  limit,  if  confirmed,  removes  the  WIMPS  proposal 
as  a  viable  option  in  solving  the  solar  neutrino  paradox. 

The  study  by  Gu,  Hill  and  Rosenwald  (1988)  of  the  asymptotic  properties 
of  the  gravity  mode  spectrum  also  found  a  fine  structure  between  asymptotic 
theory  eigenfrequencies  and  the  observed  eigenf requenci es  reported  by  Hill 
and  Gu  (1990).  This  fine  structure  has  been  interpreted  as  the  result  of 
conditions  not  being  met  for  the  applicability  of  asymptotic  theory  at  one 
or  more  radii  in  the  solar  interior.  From  an  inversion  of  the  observed  fine 
structure  by  Hill,  Gao  and  Rosenwald  (1988),  reasonably  good  agreement  is 
obtained  between  observation  and  theory  for  a  localized  perturbation  in 
internal  structure  at  r/R  =  0.23.  The  amplitude  of  the  decrease  in  the  mean 
molecular  weight  at  r/R  =  0.23  required  to  produce  the  fine  structure  has 
been  inferred  to  be  -0.9%.  Such  a  localized  decrease  in  the  mean  molecular 
weight  could  be  caused  by  a  partial  mixing  of  the  Sun's  core  in  the  past. 
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Intermediate-degree  g-modes  (those  with  angular  order  l  =  30)  were 
first  observed  in  the  late  1970's  by  Hill  and  Caudell  (1979).  However,  it 
wasn't  until  1986  that  a  preliminary  survey  was  made  of  the  1979 
differential  radius  observations  (see  3os  1982)  and  a  set  of  4  multiplets 
exhibiting  mode-locking  were  classified  by  Hill  (1986).  These  multiplets 
with  angular  order  i  *  30  and  eigenf requencies  of  =  350  pHz  were  used  as  a 
starting  point  for  the  comprehensive  analysis  discussed  in  this  work.  This 
comprehensive  study  culminated  in  the  classification  of  a  set  of  20 
intermediate-degree  g-mode  multiplets  containing  over  600  normal  modes  of 
oscillation  (Rabaey  and  Hill  1989,  1990;  Rabaey  1989a, b).  Each  of  these 
multiplets  was  found  to  contain  mode  coupled  sections.  Of  more  importance 
however  ,  is  the  internal  properties  of  the  Sun  that  can  be  inferred  from 
chis  lartte  body  of  classified  modes.  In  this  work  we  discuss  two 
significant  consequences. 

Because  these  modes  of  oscillation  are  localized  within  the  inner  50% 
of  the  Sun  by  radius  and  because  of  their  large  temperature  eigenfunctions 
implied  by  the  observed  mode-locking,  these  modes  of  oscillation  provide  a 
modification  of  the  effective  temperature  profile  defined  for  a  given 
process  in  the  Sun.  One  of  these  processes  is  the  °B  neutrino  production. 
The  second  consequence  of  these  observations  is  a  predicted  periodic 
modulation  of  the  neutrino  production  rates.  The  existence  of  a  large  set 
of  mode-coupled  gravity  modes  will  lead  to  a  low-frequency  modulation  of 
neutrino  production  rates  which  may  account  for  the  observed  periodicity  in 
the  9 B  neutrino  production  (see  Haubold  and  Gerth  1985).  The  prediction  of 
this  periodicity  in  the  neutrino  production  rates  is  unique  among  all  the 
competing  theories  that  resolve  the  Solar  neutrino  paradox  and  is  testable 
by  the  new  generation  of  Solar  neutrino  detectors. 

Theoretical  work  was  completed  which  explored  some  of  the  nonlinear 
properties  of  gravity-mode  oscillations  in  the  solar  interior  (Czarnowski 
1989).  This  work  in  particular  was  concerned  with  the  interactions  amongst 
the  nonradial  modes  of  a  single  rotationally  split  multiplet.  Motivation 
for  this  analysis  arose  both  from  general  theoretical  interest  and  from 
observations  of  solar  g-mode  oscillations,  which  show  strong  evidence  of 
nonlinear  effects  (Rabaey  and  Hill  1989,  1990;  Rabaey  1989a, b).  This  work 
provides  a  general  theoretical  framework  for  investigating  many  problems  of 
second-  and  third-order  mode  coupling  in  stellar  systems.  A  multiple-scale 
perturbation  technique  was  used.  The  formalism  presents  an  alternative  to 
that  of  Dziembowski  (1982)  and  is  more  generally  applicable. 

The  theory  provides  a  means  to  infer  the  core  amplitudes  of  the  g-modes 
exhibiting  the  nonlinear  behavior  found  by  Rabaey  and  Hill  (1989,1990)  and 
Rabaey  (1989,1990).  These  modes  have  radial  order  =  15,  angular  degree  = 
30.  This  is  a  significant  accomplishment  because  it  bypasses  the 
traditional  extrapolation  of  measured  surface  amplitudes  into  the  interior, 
which  is  a  questionable  procedure,  in  part  because  of  the  uncertainty 
regarding  the  surface  boundary  conditions  in  the  Sun.  It  was  found  that  the 
relative  radial  displacement  of  these  oscillations,  6r/r,  has  a  typical 
maximum  amplitude  in  the  interior  of  (6x10  As  a  byproduct  of  the 
analysis,  certain  characteristics  of  the  data  used  in  conjunction  with  the 
theory  allow  the  possibility  of  determining  the  linear  damping  time  of  these 
modes.  Very  preliminary  results  seem  to  indicate  a  typical  damping  time  of 
about  40  years . 
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Two  new  analysis  methods  for  solar  seismology  have  been  developed. 
Their  use  enables  both  the  solar  oscillation  frequencies  and  the  frequency 
dependencies  on  solar  internal  structure  to  be  computed  more  accurately  and 
efficiently  than  with  existing  methods.  These  two  methods  are  called  the 
continuous  orthonormalization  and  adjoint  sensitivity  methods.  There  will 
be  many  valuable  applications  of  these  two  methods  in  other  fields  of 
science  as  well  as  in  astrophysics.  This  work  was  reported  at  the  172nd 
Meeting  of  the  American  Astronomical  Society  held  in  Kansas  City,  MO,  June 
5-8,  1988  (Rosenwald  1988),  and  an  extensive  analysis  appears  in  the 
dissertation  of  Rosenwald  (1989). 

With  the  new  analysis  methods  and  the  relatively  large  number  of 
gravity  modes  classified  at  SCLERA,  the  foundation  has  been  laid  for  a  very 
effective  program  to  further  extend  the  domain  of  classified  gravity  modes. 
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of  the  AAS,  Ann  Arbor,  MI. 

June  1989  Dr.  Hill  gave  a  colloquium  at  IZMIRAN  and  Moscow 

State  University  in  Moscow,  USSR. 
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